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We apply improved nucleon-nucleon potentials up to fifth order in chiral effective field theory, along 
with a new analysis of the theoretical truncation errors, to study nucleon-deuteron (Nd) scattering 
and selected low-energy observables in ®H, ^He, and ®Li. Calculations beyond second order differ 
from experiment well outside the range of quantified uncertainties, providing truly unambiguous 
evidence for missing three-nucleon forces within the employed framework. The sizes of the required 
three-nucleon force contributions agree well with expectations based on Weinberg’s power counting. 
We identify the energy range in elastic Nd scattering best suited to study three-nucleon force effects 
and estimate the achievable accuracy of theoretical predictions for various observables. 
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Chiral effective field theory (EFT) provides a pow¬ 
erful framework for analyzing low-energy nuclear struc¬ 
ture and reactions in harmony with the symmetries of 
quantum chromodynamics (QCD), the underlying the¬ 
ory of the strong interactions. It allows one to derive 
nuclear forces and currents in a systematically improv¬ 
able way in terms of a perturbative expansion in powers 
of Q S (p/Ab, M,r/At,), the so-called chiral expansion. 
Here, p refers to the magnitude of the nucleon three- 
momentum, is the pion mass and Ab is the break¬ 
down scale of chiral EFT [T]. One finds, in particular, 
that the leading-order (LO) contribution to the Hamil¬ 
tonian at order Q® and the first corrections at order Q® 
(NLO) are given solely by nucleon-nucleon (NN) oper¬ 
ators while three-nucleon forces (3NFs) appear first at 
order Q® (N®LO) (see [2] and references therein). Four- 
nucleon forces are even more suppressed and start con¬ 
tributing at order Q'^ (N®LO). The chiral power counting 
thus provides a natural explanation of the observed hier¬ 
archy of nuclear forces. 

With accurate N®LO NN potentials being available 
since about a decade 011], the main focus of research 
has moved in recent years towards the 3NF 00. While 
providing a small corrections to the nuclear Hamiltonian 
as compared to the dominant NN force, its inclusion 
is mandatory for quantitative understanding of nuclear 


structure and reactions. Historically, the importance of 
the 3NF has been pointed out already in the thirties 0 
while the first phenomenological 3NF models date back 
to the fifties. However, in spite of extensive efforts, the 
spin structure of the 3NF is still poorly understood 0 . 

Chiral EFT is expected to provide a suitable theoreti¬ 
cal resolution to the long-standing 3NF problem. Indeed, 
the leading chiral 3NF has already been extensively ex¬ 
plored in ah initio calculations by various groups and 
found to yield promising results for nuclear structure and 
reactions 00- The first corrections to the 3NF at order 
(N®LO) have also been derived [5HTT] (and appear 
to be parameter-free) while the sub-subleading contribu¬ 
tions at order Q® (N'^LO) are being derived 

On the other hand, understanding and validating the 
fine details of the 3NF clearly requires precise and sys¬ 
tematic NN potentials and a reliable approach for esti¬ 
mating the accuracy of theoretical predictions at a given 
chiral order. Refs. [m [10 initiated the direction that 
we follow here by developing a new generation of chi¬ 
ral EFT NN forces up to N'^LO, in which the amount 
of finite-cutoff artefacts has been substantially reduced 
by employing a novel ultraviolet regularization scheme, 
and by introducing a new procedure for estimating the 
theoretical uncertainty. 

In this Letter we, for the first time, apply these novel 
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chiral NN forces beyond the two-nucleon system and 
demonstrate their suitability for modern ab initio few- 
and many-body methods. By applying the new method 
for error analysis, we present unambiguous evidence for 
missing 3NF effects and demonstrate that the size of the 
required 3NF contributions agrees well with expectations 
based on Weinberg’s power counting. We also estimate 
the theoretical accuracy for various observables achiev¬ 
able at N^LO and identify the energy region in elastic 
Nd scattering that is best suited for testing the chiral 
3NF. 

We first describe our procedure for estimating the the¬ 
oretical uncertainty. Let X{p) be some observable with 
p referring to the corresponding momentum scale and 
i = 0,2,3,..., a prediction at order Q* in the 
chiral expansion. We further define the order-Q* correc¬ 
tions to X{p) via 

AxW = X^^'>-X^^-^\ i>3, (1) 

so that the chiral expansion for X takes the form 

x(^) = -p AX(2) -p ... -p AA(*) . (2) 

Generally, the size of the corrections is expected to be 

AAW = (3) 

In |16] . the validity of this estimation was confirmed 
for the total neutron-proton cross section. In Refs. 
[13 [S], quantitative estimates of the theoretical uncer¬ 
tainty dXh) of the chiral EFT prediction were made 
using the expected and actual sizes of higher-order con¬ 
tributions. Specifically, the following procedure was em¬ 
ployed: 

(4) 

= max(Q*+i|A(°)|, \AX^^^\y 2<j<i 

where i>2 and Q = max(p/A;,, M,r/A{,) with Af, = 600, 
500 and 400 MeV for the regulator choices of i? = 
0.8 — 1.0 fm, R = 1.1 fm and R = 1.2 fm, respec¬ 
tively. The sizes of actual higher-order calculations pro¬ 
vide additional information on the theoretical uncertain¬ 
ties, which we use by adding the conditions 

> max (5) 

to estimates of lower-order uncertainties. 

The above procedure for estimating the uncertainty 
needs to be adjusted in order to account for the neglect 
of many-body forces in the present analysis. In partic¬ 
ular, iterating the NN T-matrix in the Faddeev equa¬ 
tion generates contributions whose short-range behav¬ 
ior is order- and regulator-dependent. For low-energy 
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FIG. 1: (Color online) Chiral expansion of the EgB based 
on the NN potentials of Refs. [13 US] for the regulator R = 
1.0 fm and using Q = M-n/X- Left (middle) panel shows 
incomplete results based on NN forces only, with uncertainties 
being estimated via Eqs. (|^[^ (Eqs. ([^[^ for i > 3). Right 
panel shows the projected results assuming that the LECs in 
the N^LO 3NF are tuned to reproduce the Eg^ and using 
Eqs. (|^[^ to specify the uncertainty. 

Nd observables calculated in the EFT framework, ap¬ 
proximate scheme independence is restored upon per¬ 
forming renormalization, i.e. upon expressing the bare 
low-energy constants (LECs) accompanying short-range 
3NEs at orders Q®, ... in terms of observable quan¬ 
tities, such as the triton binding energy. In practice, this 
is achieved by htting the corresponding LECs to exper¬ 
imental data. Therefore, when performing incomplete 
calculations based on NN interactions only, the estima¬ 
tion in Eq. ([^ is not justified at or beyond N^LO, the 
chiral order at which the contact 3NF starts contributing. 
We, therefore, adopt here a slightly modihed procedure 
for estimating the uncertainty for i > 3, namely 

(5A(*) =max(Q*+i|A(°)|, Q^-'^\AX<~‘^'>\, Q*-2|AA(3)|^, 

( 6 ) 

and do not employ Eq. (§. However, to be conservative 
in our estimates, we further require that 

(5A(2) > ^ (7) 

The dependence of the chiral NN forces on the local 
regulator R over the range 0.8 ... 1.2 fm has been exten¬ 
sively investigated in Ref. [13 showing that cutoff arti¬ 
facts become visible for R > 1.0 fm. On the other hand, 
we seek to obtain many-body results as close to conver¬ 
gence as possible, and this favors the largest feasible value 
of R. We therefore balance these competing conditions 
with the choice of i? = 1.0 fm in this work. 

Our results for the chiral expansion of the ground 
state energy (Egg) using Q = are visualized in 

Fig.m Assuming that the LECs which enter the short- 
range part of the 3NE can be tuned to reproduce the 
Egg, we can already at this stage present a complete result 
up to N'^LO for the chiral expansion of this observable, 
see the right panel in Fig. [^ As expected, we observe 
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FIG. 2: (Color online) Predictions for Nd total cross section 
based on the NN potentials of Refs. |15[ I16| for R = 1.0 fm 
without including the 3NF. Theoretical uncertainties (blue) 
are estimated via Eqs. 0 and 0 for chiral order i = 0, 2 and 
via Eqs. 0 and Q for i > 3. Experimental data are from 
Ref. [T7| . 

that Eqs. 00 provide a more reliable approach for er¬ 
ror estimation in calculations based on NN interactions 
only, while using Eqs. 0 amounts to overestimating 
the actual error. The I^LO (N'^LO) results for the 
Egs are expected to be accurate at the level of ^50 keV 
(^10 keV) for the regulator choices of i? = 0.8, 0.9 and 
1.0 fm. Note that the size of the 3NE contribution agrees 
well with the uncertainty at NLO, which reflects the es¬ 
timated impact of the N^LO contributions to the Hamil¬ 
tonian. This is fully in line with expectations based on 
the Weinberg power counting mm- 

We now turn to Nd scattering observables. Our predic¬ 
tions for the Nd total cross section are visualized in Eig.[^ 
Similar to the Egg, one observes a significant discrep¬ 
ancy between the theoretical predictions based on the 
NN forces only and data, which provides clear evidence 
for missing 3NF contributions. The size of the discrep¬ 
ancy agrees within 1.5 times the estimated size of N^LO 
corrections shown by the NLO error bars. Interestingly, 
the discrepancy at the lowest energy of 10 MeV is much 
smaller than the estimated size of N^LO contributions. 
Given that the cross section at low energy is governed 
by the S-wave spin-doublet and spin-quartet Nd scatter¬ 
ing lengths, this observation can be naturally explained. 
Indeed, the spin-quartet scattering length is almost an 
order of magnitude larger than that of the spin-doublet 
and much less sensitive to the 3NF as a consequence of 
the Pauli principle. 

Our predictions for Nd differential cross section and an¬ 
alyzing powers Ay{N), Ayy and A^^ are shown in Figs.0 
1^ At the lowest energy of 10 MeV, there is little appar¬ 
ent need for 3NF effects except for Ay. Interestingly, the 
fine-tuning nature of this observable is clearly reflected in 
large theoretical uncertainties at NLO and N^LO. Start¬ 
ing from En = 70 MeV, one observes clear discrepancies 
between our predictions and data for the cross section 
and tensor analyzing powers which are expected to be 
explained by the 3NF. In all cases, the required 3NF 
contributions are of a natural size. Based on the width 





FIG. 3: (Golor online) Predictions for the differential cross 
section and nucleon Ay in elastic Nd scattering based on the 
NN potentials of Refs. |151116| for R = 1.0 fm without in¬ 
cluding the 3NF. Theoretical uncertainties are estimated via 
Eqs. 0 and 0 for chiral order i = 2 and via Eqs. 0 and 
0 for i > 3. The bands of increasing width show estimated 
theoretical uncertainty at N^LO (red), N^LO (blue), N^LO 
(green) and NLO (yellow). The dotted (dashed) lines show 
the results based on the GD Bonn NN potential [18] (CD Bonn 
NN potential in combination with the Tucson-Melbourne 3NF 
[19jL For references to proton-nucleon data (symbols) see 
Ref. 1^. 

of the bands, one may expect Nd scattering observables 
at N^LO to be accurately described up to energies of 
at least 200 MeV. It is also comforting to see that the 
accuracy of chiral EFT predictions for Nd and NN [TS] 
scattering observables at the same energy is comparable. 
We further emphasize that the improved NN potentials 
of Refs. [iSKIll show clearly smaller finite-cutoff artifacts 
as compared to the N^LO potentials of Refs. [31 Hj and, in 
particular, do not lead to distortions in the cross section 
minimum that were found in Ref. |20j . 

Next, we apply the improved NN potentials to A > 3 
systems. We present in Fig. order-by-order calcula¬ 
tions of selected observables for ^He and ®Li. The re- 
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FIG. 4: (Color online) Predictions for the tensor analyzing 
powers Ayy and in elastic Nd scattering based on the NN 
potentials of Refs. [13 US] for i? = 1.0 fm without including 
the 3NF. For notations see Fig. [^ 

suits for ^He are obtained both by solving the Faddeev- 
Yakubovsky (FY) equations and with the no-core shell 
model (NCSM) |8], which agree to within the estimated 
uncertainties of these methods. The numerical uncer¬ 
tainties in the FY solutions are a few keV for the energy 
and about 0.001 fm for the point-proton radius (cp). The 
numerical uncertainties from incomplete convergence of 
the NCSM (see Ref. [21] for details) are shown as error 
bars (color online: red) together with the estimated theo¬ 
retical uncertainties from the truncated chiral expansion 
with Q = Mtt/A f, (color online: blue). 

For the ®Li energies, we carried out Similarity Renor¬ 
malization Group (SRG) evolution [^ in order to en¬ 
hance the convergence rate of the NCSM calculations 
that were performed in basis spaces up through A^max = 
12 and extrapolated to the infinite matrix limit following 
Ref. [22]. We retained the induced 3NF arising from the 
SRG evolution, see Ref. [53] for details, and this produces 
results for the ®Li energies in Fig.j^that are independent 
of the SRG scale over the range a = 0.04 — 0.08 fm^ to 
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FIG. 5: (Color online) Predictions for EgB and Cp of '’He and 
the energies of the lowest two states of ®Li based on the NN 
potentials of Refs. [13 [13 for R — 1.0 fm without including 
the 3NF. Theoretical uncertainties (blue) are estimated via 
Eqs. and <© for chiral order i = 0, 2 and via Eqs. <§ and 
0 for i > 3. Numerical uncertainties from the NCSM (red) 
are estimated following Ref. m- 

within our quoted many-body uncertainties. For exam¬ 
ple, at N'’L0 we obtain Egg = —26.9(4) (—26.9(2)) MeV 
at a = 0.04(0.08) fm^ for ®Li where the quantified numer¬ 
ical uncertainty in the last digit of the energy is quoted 
in parenthesis. 

The patterns for the energies in Fig. as well as for 
the Tp of ^He are very similar to the pattern for the Egg 
of in Fig. and the Nd total cross section at 10 MeV 
in Fig. As in ^H, we again observe underbinding in¬ 
dicative of the need for 3NFs, especially at N^LO and 
N^LO. This underbinding is correlated with larger rp in 
^He, which is expected to decrease toward the experi¬ 
mental result as Egg is lowered toward experiment with 
the inclusion of 3NFs. Note that the energy of the first 
excited state in ®Li, with = d’*', follows the same pat¬ 
tern as the ground state energy, leading to an excitation 
energy that depends much less on the chiral order than 
one might naively expect based on the theoretical uncer¬ 
tainties of the binding energies. 

To summarize, we have studied in this Letter selected 
few-nucleon observables using improved chiral NN po¬ 
tentials of Refs. [Iil[I3 up to N'’L0. Our results suggest 
that these new chiral forces are well suited for modern 
ab initio few- and many-body methods. Using the novel 
approach for error analysis introduced in Ref. |15j . we 
found truly unambiguous evidence for missing 3NF ef¬ 
fects by observing discrepancies between our predictions 
and experimental data well outside the range of quanti¬ 
fied uncertainties. The magnitude of these discrepancies 
is found to match well with the expected size of the chi¬ 
ral 3NF whose dominant contribution appears at N^LO. 
Furthermore, we have demonstrated that the predictions 
for Nd and NN scattering observables at the same energy 
have comparable accuracy, in agreement with the general 
principles of EFT. Most importantly, the expected theo¬ 
retical uncertainty for Nd scattering observables at N^LO 
and N'’L0 in the energy range of Eiab — 70 — 200 MeV is 
shown to be substantially smaller than the observed dis- 
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crepancies between state-of-the-art calculations and ex¬ 
perimental data. This suggests that chiral EFT at these 
orders should be capable of resolving the long standing 
3NF problem in nuclear physics. Work on the explicit 
inclusion of the consistent 3NFs is in progress. 
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